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Polymer chains adsorbed onto oppositely charged spherical colloidal particles can significantly 
modify the particle-particle interactions. For sufficient amounts of added polymers, the original 
electrostatic repulsion can even turn into an effective attraction and relatively large kinetically stable 
aggregates can form. The attractive interaction contribution between two oppositely particles arises 
from the correlated adsorption of polyions at the oppositely charged particle surfaces, resulting in 
a non-homogeneous surface charge distribution. Here, we investigate the aggregation kinetics of 
polyion-induced colloidal complexes through Monte Carlo simulation, in which the effect of charge 
anisotropy is taken into account by a DLVO-like intra-particle potential, as recentely proposed by 
Velegol and Thwar [D, Velegol and P.K. Thwar, Langmuir, 17, 2001]. The results reveal that in the 
presence of a charge heterogeneity the aggregation process slows down due to the progressive increase 
of the potential barrier height upon clustering. Within this framework, the experimentally observed 
cluster phases in polyelectrolyte-liposomes solutions should be considered as a kinetic arrested state. 
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1. INTRODUCTION 

The addition of oppositely charged polyions to a sus- 
pension of charged colloidal particles gives rise to an in- 
triguing and partially unexpected phenoinenologyii2,j re- 
sulting in the formation of aggregates which play an im- 
portant role in a wide range of implications^i^i^, from 



&J.&. 



to 



membrane biophysics and soft matter physics^ii 
biotechnological processes, such as therapeutic delivery 
systems^. 

These aggregates are governed by a delicate balance 
between attractive and repulsive interactions resulting in 
the appearance of an equilibrium (or a kinetically ar- 
rested) cluster phase, where single particles, stuck to- 
gether by an electrostatic "glue" acting between oppo- 
sitely charge domains, form reversible, relatively large, 
complexes. In this system, the short-range attraction 
contribution is promoted by the addition of " adsorbing" 
polyions which form a two-dimensional strongly corre- 
lated, short-range ordered, structure on the surface of 
the oppositely charged particle, contrarily to what hap- 
pens for " non- adsorbing" polymers, where the attraction 
contribution is produced by the unbalanced osmotic pres- 
sure in the depletion regime. 

In fact, correlated adsorption of polyion chains on the 
surface of oppositely charged particles induces two dif- 
ferent phenomena: "charge inversion"i^i^ii and "reen- 
trant condensation" ^i^i^° . The first one occurs when a 
colloidal particle binds several oppositely charged multi- 
valent ions (polyions), so that its effective charge inverts 
its sign. The second effect, concomitant to the charge 
inversion, consists in the formation of particle aggregates 
whose average size increases on increasing the polyion 
concentration, until it reaches a maximum (at the point 
of charge inversion), decreasing afterwards to the initial 
value. 

This phenomenology has been investigated in a vari- 
ety of different colloidal systems both for particles bear- 



ing a negative charge (such as DOTAP liposomes) in the 
presence of anionic polyions^ and for positively charged 
particles (such as hybrid niosomcs) in the presence of 
cationic polyions^ ^. The occurrence of a sharp reentrant 
condensation on increasing the adsorption of polyion on 
charged liposomes has been documented in a series of 
article a^'^'^" . It has been shown that the size of the 
aggregates goes through a maximum and their over- 
all charge inverts its sign in concomitance to the point 
of charge inversion, as briefly reviewed in Fig. [1] for 
some typical cationic or anionic colloidal particle sys- 
tems. The formation of aggregates and, consequently, 
the stability of colloidal charged particles (solid parti- 
cles or soft vesicles) is usually satisfactorily explained 
in the framework of Derjaguin-Landau-Verwey-Overbeek 
(DLVO) theoryiS,, which accounts for electrostatic and 
van der Waals interactions between two approaching par- 
ticles. In the traditional DLVO theory, the particle sur- 
face is assumed to be uniformly charged. In the case 
of polyion-correlated adsorption onto the surface of the 
particle, this assumption is no longer valid and the usual 
DLVO theory might not account for the effects arising 
from the surface charge heterogeneity. Indeed, the ad- 
sorption of the polymer generates an ordered charged 
structure, strongly dependent on the polyion size and on 
the valence and linear charge density^ . Some experimen- 
tal evidences of such a structure can be found in Refs.[13] 
and [14]. Transmission electron microscopy [TEM] mea- 
surements suggest the formation of a local "patchy charge 
distribution". In this respect, it is reasonable to expect 
that attractive interactions can arise when "counterion 
domains" on one particle align to a " counterion- free do- 
mains" on another approaching particle. 

Recently, Velegol and Thwar^^ have developed an an- 
alytical model for the effective particle-particle potential 
($) in the case of non-uniform charge distributions, based 
on the Derjaguin approximation and on an extention of 
the Hogg-Healy-Fuerstenau [HHF]^^ model for randomly 
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FIG. 1: Fig. 1 - A brief review of the experimental evi- 
dence showing a non-monotonic increase of the size of the 
polyion-induced charged particle aggregates, in the presence 
of oppositely charged polyions. Both positively and nega- 
tively charged particles are shown. Left panels: reentrant 
condensation (aggregate radius R) (A) and charge inversion 
(^-potential) (B) of positively charged liposomes, in the pres- 
ence of negatively charged polyions. Right panels: reentrant 
condensation (aggregate size R) (A) and charge inversion 
(("-potential) (B) of negatively charged hybrid niosomes, in 
the presence of positively charged polyions. The data are 
shown as a function of the charge molar ratio ^ normalized 
to the value ^o at which the measured ("-potential goes to 
zero. ^ is defined as the polyion to the lipid molar ratio. 
The experimental values (hydrodynamic radius measured by 
means of dynamic light scattering and (-potential measured 
by means of Doppler laser electrophoretic techniques) refer to 
different colloidal systems. Positive charged particles: (v)- 
DOTAP liposomes (0.8 mg/ml) and polyacrylate sodium salt; 
(•): DOTAP liposomes (1.7 mg/ml) and polyacrylate sodium 
sallii^. Negative charge particles: (o): hybrid niosomes and 
a-polylysine; (A): hybrid niosomes and e-polylysine; (D): hy- 
brid niosomes and PE'VP (ionization degree 65%); (o): hybrid 
niosomes and PEVP (ionization degree 95%)^^. Hybrid nio- 
somes are built up by Tween20, Cholesterol and dicethylphos- 
phate and the cationic polyion PEVP is Poly [N-ethy 1-4- vinyl 
pyridinium] bromide. 



charged surfaces. The resulting spherically symmetric 
potential depends on the values of the average surface 
potential (here assumed equal to the (^-potential) and on 
the standard deviation a of the surface potential among 
different regions on the particle surface. In particular, 
Velegol and Thwar^^ showed that the value of the second 
moment of the charge distribution lowers the electrostatic 
repulsion. 

In this work, by means of MC simulation, we ana- 
lyze the aggregation process of spherical charged particles 
interacting via the Velegol and Thwar^^ DLVO-like po- 
tential. Starting from isolated monomers, we follow the 
growth of the aggregates until a kinetically arrested state 



is reached. We study the final size of the aggregates for 
several values of the ("-potential and its standard devia- 
tion a. The results of the simulations confirm that the 
final size of the aggregates grows on increasing a and 
on decreasing (", providing a microscopic interpretation 
of the re-entrant condensation phenomenon of polyion- 
induced charged colloidal particle aggregation (Fig.[T|). 



2. THEORETICAL BACKGROUND 

According to the HHF model, the electrostatic po- 
tential of mean force 4> between two dielectric colloidal 
spherical particles of radius Ra and Rb and surface po- 
tentials ipA and ipB , respectively, is given by 
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€itRaRb 
Ra + Rb 



[(^ 



^-rV^B)ln(l-e" 



-2kH\ 



+ 2^AipB In (coth ^)] 



(1) 



Here, e is the dielectric permittivity of dispersing 
medium, k~^ is the Debye screening length and H is the 
minimal distance between the two particle surfaces (see 

fig-H- 

Eq. ^ , which is derived under the conditions kR^ S> 1 




FIG. 2: "Variables involved in the Derjaguin approximation, 
for two approaching particles of radius Ra and Rb, respec- 
tively 



and H/Ra ^ 1 (a=A,B), holds under the following as- 
sumption: i) ipA and ipB are relatively small (less than 
25 mV at room temperature); ii) the particles share the 
same chemical nature and are charged by similar mech- 
anism; iii) the Derjaguin approximation holds^^. 

Velegol and Thwar— assume that the surface of a par- 
ticle a is partitioned in N regions (labelled by i, with 
i = 1, ..., Na) of area 5", each of them characterized by a 
different value of the surface potential -0^ . These regions 
must be of sufficient size L ~ y/S, so that lateral interac- 
tions can be considered negligible (practically this means 
that L ^ K^^ and L <C H). The surface potentials ■0A 
and ipB of two approaching particles in the HHF theory 
are thus replaced by the random value of the potential 
-0^ and 0g, where i and j label the two domains (on 
opposite particles) facing each other. Defining ("„ as the 



average potential, one can write ■i/)^ = Ca + ^^ where S^^ 
a random contribution that varies among the different 
regions. The (5^ are independently distributed and not 
correlated, i.e. 



(CV-^) = 



CaCb, 



a+' 



'.S'\ 



ii a = [3 = A 
HA^B 
iia = (3 = B 



(2) 



where 6^^ is the Kroenecker delta function and aa the 
standard deviation of surface potential on spheres of type 
a. 

The expression for the resulting potential of mean force 
provided by Velegol and Thwar— is 

M 

•^ - YJ2i(^A+Cl+2CASf+2CBSf+S\+Sl){l-cothKh,) 



H 

s 

e 

V 



4- 

2- 







•') .■ ^ ^.-"^^ ■ — - — ^ 


2- 
4- 
6- 


f.!/ 

I'/ 
'■■'/ 

iL-^ \ , \ 1 \ 1 \ 1 1 



10 15 

H [nm] 



20 



25 



2{CA + 6f){CB + Sn 
cosh nhj 



}S 



(3) 



where the sum runs over the M facing regions and hi 
is the gap between jth regions. The ensemble averaging 
gives: 

(4) 
In order to evaluate the sum in Eq. ^, the Derjaguin 
approximation is applied (Fig. [J), i.e. 



, ^ Ti , Ra + Rb 2 
2RaRb 



(5) 



and, for any function F{h), the following expression are 
employed 



AI 



Y,Fih,)A,^2^ 
1=1 



RaRb 
Ra + Rb j 



F{h)dh (6) 



where r is the flat ring radius (see Fig. [2). 

The resulting mean force pair interaction potential is 

given by 

($> - l^fll [(Cl + Cl + <r\ + 4) ln(l - e--^)+ 



FIG. 3: Mean force potential profiles between two identical 
spheres calculated from Eq. ([T]) for different values of C^a = 
C,B = C (from 11 to 19 mV) and for a constant value of a a = 
(7_B = (J=15 mV. Curves are plotted in units of the thermal 
energy ki,T at room temperature for Ra = Rb = R — 40 nm 
and K~^ = 10 nm. Solid line: (" = 11 mV; dashed line: (" = 13 
mV; dotted line: (^ = 15 mV; dot-dashed line: ^ = 17 mV; 
dot-dot-dashed line: (^ = 19 mV. 
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The mean force potential for typical values of the (^- 
potential and standard deviation a is shown in Figs. [3] 
and m The potential combines a repulsive net charge- 
dependent monopole term (Cq ^ 0) and an attractive 
multipole term (ctq, 7^ 0) arising from the presence of 
random charge heterogeneity on the particle surface. 



FIG. 4: Mean force potential profiles between two identical 
spheres calculated from Eq. (0 for different values of aA ~ 
Ob ~ cr (from 5 to 25 mV) and for a constant value of C^a ~ 
C^B ~ C=15 mV. Curves are plotted in units of the thermal 
energy kbT at room temperature for Ra = Rb = R = 40 
nm and k~^ = 10 nm. Solid line: a = 25 mV; dashed line: 
(T = 20 mV; dotted line: cr = 15 mV; dot-dashed line: it = 10 
mV; dot-dot-dashed line: a — 5 mV. 



These two components lead up to a global maximum 
and to the presence of a potential barrier. In the case 
of particles with the same net charge, an attractive com- 
ponent is always present since the first term in Eq. ([7]) 
is negative. If the two identical particles are uniformly 
charged (aa — 0), both the attractive component and the 
global maximum vanish, recovering the HHF expression 
for two identical spheres (Eq. [1]). 

The height of the potential barrier that two approach- 
ing particles must overcome in order to aggregate and 
the distance between the particle surface at which this 
maximum occurs can be evaluated from Eq. ([7]). For 
two identical particles {Ra = Rb = R) we obtain: 
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We point out that: i) the height of potential barrier 
does not depend on Debye screening length k~^. Con- 
sequently, addition of simple salt in bulk phase (varying 
the ionic strength of the solvent) does not change the 
strength of the interaction but only modifies the distance 
at which this maximum interaction occurs; ii) using Der- 
jaguin approximation, a linear dependence of the barrier 
height on the radius R of the particles arises, underlying 
that larger particles interact more intensely than smaller 
ones and thus aggregation in the presence of large clus- 
ters is relatively inhibited; in) ($)ma£c is zero when C = 
and limo.^o+ Hmax = 0, that is consistent with the state- 
ment that there is no attractive component for uniformly 
charged particles and that a lower value of a determines 
a short-ranged attractive interaction. 

In the framework of the above stated model, the ag- 
gregation process can be viewed as a thermally activated 
process strictly linked to electrostatic parameters (the C,- 
potential and its standard deviation a) and to the clus- 
ters dimensions. Colloidal particles with smaller curva- 
ture radius, lower ^-potential or higher charge anisotropy 
will be characterized by a faster aggregation dynamics. 
The aggregation process will slow down due to the pro- 
gressive increase of the repulsive barrier on increasing the 
size of the aggregates. 

In the next section, we introduce a MC simulation 
method to analyze the aggregation dynamics of spher- 
ical clusters and to evaluate the evolution of the mean 
cluster radius for several values of C-potential and stan- 
dard deviation a, chosen to model the corresponding ex- 
perimental range of values observed for polyion-induced 
liposome cluster aggregation. 



In our calculations, the contribution to the attraction 
associated to the van der Waals interaction between two 
approaching particles is neglected. The reason why such 
interaction can be safely neglected arises from the typ- 
ical range associated to the van der Waals interaction, 
which dyes off before Hmax- Indeed, for liposome par- 
ticles (aqueous cores of radii i?i and i?2, respectively, 
covered by phospholipid shell of thickness d) , the van der 
Waals interaction Vydw can be written a&ii 



VvdW — 
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H + 2d H + d 



H{H + 2d) 



1 
H 



2{H + d) 



(10) 



(where A is the Hamaker constant). For typical values 
of the parameters involved, eq. [10] results in a range of 
attraction which becomes negligible beyond 1 -;- 2 nm. As 
can be seen from Figs. [3] and [H this distance is always 
smaller than Hmax ■ Consequently, we neglect the van der 
Waals interaction in the irreversible aggregation process. 



3. SIMULATION 

We study a system composed by Np = 10000 colloidal 
particles of initial diameter 2R = 80 nm in a cubic box 
of volume V with packing fraction = AnpR^ /'i = 0.01 
where p = Np/V is the number density. We carry out 
MC simulation using local metropolis algorithm at T = 
298 K. Particles interact via a short range potential 
defined by Eq. ([7]). We study five different pairs of ^ 
and a values, comparable to typical values measured in 
liposome solution£^*''ii. We choose 1/k = 10 nm. We 
recall that the value of screening length does not modify 
the barrier height. 

Because of the compactness of the formed clusters, 
aggregation is modeled as an "oil drop-like process" in 
which the shape of the aggregates retains a spherical 
form. More precisely, when two approaching particles 
{A and B) overcome the potential barrier (surface gap 
distance H < Hmax), they aggregate forming an unique 

particle with radius R = (^R\ + Rg) and positioned 
in the center of mass of two aggregating spheres. We also 
assume that S (the size of the uniform potential regions 
on the particle surface in the Velegol-Thwar expression) 
is independent of the cluster size. To incorporate a Brow- 
nian dynamics, in the MC algorithm, the ith particle is 
selected with a probability proportional to Rq/ Ri, where 
Rq is the initial radius and Ri is the radius of the i- 
th aggregatei^iiSiSS. Each selected aggregate is moved 
in each direction by a random quantity uniformly dis- 
tributed between ±0.2 nm. Simulations have been car- 
ried out by varying both the ^-potential and the variance 
a^ . The ^-potential values have been chosen similar to 
the typical ones measured in different colloidal systems 
(see Fig. 1 and Refs. 10-11). On the contrary, values for 



(T^, in absence of any experimental indication, have been 
chosen arbitrarily, within the validity of the model and in 
any case within a range of reasonable values. The aggre- 
gation process progressively slows down and simulations 
are interrupted when a plateau in the time dependence 
of the aggregate average radius or mass is reached. 



4. RESULTS AND DISCUSSION 

In Figs. [n]and[ni we present results of the aggregation 
process investigated for different values of the character- 
istic parameters (^-potential and standard deviation a). 

First, we will discuss the time evolution of mean clus- 
ter radius at constant a on varying the ^-potential. Fig. [6] 
shows that on increasing of the C-potential, the aggrega- 
tion process slows down and that the aggregates reach 
at long time a final limiting size. The slowing down of 
the dynamics is already seen at the early stage of the 
aggregation process, as shown in the inset of Fig. [6l con- 
sistent with the effect of C-potential on the height of the 
potential barrier. 

Fig. [5] shows the complementary case, i.e, the effect of 
the different values of the standard deviation cr at a con- 
stant value of the ^-potential. Again, data confirm that 
the increase of the random fluctuations of the charges on 
the particle surface speeds up the aggregation process. 
In all the cases investigated, the growth process slows 
down at long times and the aggregates appear to reach 
a long-time limit value, providing a strong evidence for 
a dynamics slow down and an arrest on the timescale 
sampled by our simulation. This finding suggests that a 
kinetically arrested state can be generated in this class 
of systems by the coupling between the aggregates size 
and the electrostatic barrier. An estimate of the charac- 
teristic size reached by the aggregates at long time can 
be calculated assuming that the arrest is observed on the 
explored time scales when the interaction energy barrier 
reaches an ad-hoc value. Indeed, inverting Eq. [51 the re- 
lation connecting the radius of interacting aggregates and 
the potential barrier is easily obtained. We find that, in 
all the cases investigated, for all the couples of C-potential 
and standard deviation a we have considered, the plateau 
of the average aggregate size corresponds to a characteris- 
tic barrier {^)max of the order of about 10 ki,T (see Figs. 
[S] and [3 where the normalized size < R > /Rq-I calcu- 
lated from eq. [8] for different values of the C-potential 
at fixed a or conversely, for different values of a at fixed 
C-potential is shown). 

It is worth nothing that C a-nd a values are consid- 
ered to be both constant during the aggregation process 
that evolves from the single polyion-decorated liposome 
to liposome cluster. This assumption is justified as a con- 
sequence of the local character of the interaction, where 
only a local charge (potential) distribution determines 
the repulsive and attractive components of the pair po- 
tential. The distribution of the local charge of the ele- 
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FIG. 5: Some typical time evolution of normalized mean clus- 
ter radius < R > /Ro-1- Simulations have been carried out 
for different values of the standard deviation a (in the range 
from 5 to 25 mV) with a constant value of the C-potential C= 
15 mV. The inset shows, in linear scale, the evolution at short 
times of the mean size of the aggregates, for different values 
of the standard deviation a (cr=25, 20, 15, 10, 5 mV, from 
top to bottom, respectively). 



mentary unit (single liposome) surfaces can be considered 
independent of the size of clusters. 

In all the simulations, at short time below the slow 
down, the average size of clusters increases according to 
a power law {R)/Ro — 1 ~ at'' with a pre- factor that de- 
creases from a = 7.2 • 10~® to a = 1.6 • 10~^. At the same 
time, the exponent b increases from b ~ 0.70 to 0.88, as 
the C-potential varies from C = H mV to C = 19 mV at a 
constant value of cr = 15 mV. The same behavior charac- 
terizes the a dependence at fixed C = 15 mV in the range 
from a — 25 mV down to cr = 5 mV. As expected, the 
aggregation process does not follow the diffusion limited 
cluster aggregation (DLCA model^i) due to the presence 
of significant potential barriers controlling the rate of ag- 
gregation. Deviation from DLCA is also evidenced by 
the not negligible potential barrier height at which the 
arrested state occurs (of order of 10 kbT). 

Fig. [7] shows the cluster mass distributions at differ- 
ent stages, during the aggregation process for the case of 
the couple of parameters C = 15 mV and cr = 15 mV. 
The system evolves from the initial mono-disperse dis- 
tribution toward a broader distribution characterized by 
a well defined peak, whose position increases with time. 
The distribution freezes when the system kinetically ar- 
rests. The inset in Fig. [7] shows the same distribution in 
a semi-log scale to highlight the presence of a distribution 
tail compatible with an exponential distribution. Similar 
behavior is observed for all the cases investigated. The 
main differences concern the shape of the distributions 
in the arrested states, some of which shown in Fig. [S] 
(see also Table HI . The final size distribution has a 
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FIG. 6: Some typical time evolution of normalized mean clus- 
ter radius < R > /Ro-i- Simulations have been carried out 
for different values of ^-potential (in the range from 11 to 19 
mV) with a constant value of the standard deviation a= 15 
mV. The inset shows, in a linear scale, the evolution at short 
times of the mean size of the aggregates for different values 
of the (^-potential (C= 11, 13, 15, 17, 19 mV, from top to 
bottom, respectively). 
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FIG. 7: Cluster size distribution of the mass M of the ag- 
gregating particles normalized to the mass Afo of the initial 
particles at different aggregation states for (^ = 15 mV and 
cr = 15 mV. (o): (M)/Mo = 2; (A): {M)/Mo = 5; (v): 
{M)/Mo = 10; {<»: {M)/Mo = 25; (•): (M)/Mo = 50. In- 
set: the same distributions in a semi-log scale. 



smaller standard deviation aM in the case of higher val- 
ues of the C-potential or lower values of a. The average 
mass instead shows an opposite trend, i.e., it increases 
for smaller values of the ^-potential or higher values of 
a. 

A possible explanation of the influence of the electro- 
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TABLE I: Standard deviation aM of mass distributions when 
{M)/Mo = 50 for the different electrostatic parameters (("- 
potential and its standard deviation a) used in simulations. 
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FIG. 8: (Color on line) Cluster mass distribution in a log- 
scale for the cluster distributions closest to the arrested state; 
(.): C=ll mV, cr=15 mV; (a): C=15 mV, cr = 15 mV; (T): 
C=15 mV, a = 25 mV; (♦): C=15 mV, a = 20 mV; (■): C=13 
mV, cr = 15 mV . 



static parameters, ^-potential and standard deviation cr, 
on the variance cta/ of the cluster size distribution is pro- 
vided in Fig. [HI Here, we compare the potential profiles 
(eq. [7]) where a particle interacts with particles of differ- 
ent sizes, for two different values of the ^-potential (11 
and 15 mV, respectively). The results show that there 
is a preferential aggregation which favors the increase of 
smaller cluster rather than the cluster of larger size. In- 
deed, the inter-particle potential (Eq. ([7])) predicts that 
the test particle (particle of mass lOA^fo in the example 
sketched in Fig. ^ preferentially interact with clusters 
of smaller size (particle of mass 50Afo) rather than with 
clusters of larger size (particle of mass 200Mo) This kind 
of preferential interaction, which becomes relatively more 
and more favorable for larger values of (, favors the aggre- 
gation of the smaller clusters in the system. These priv- 
ileged interactions result in the formation of aggregate 
size distributions whose width decreases with increasing 
(^ or decreasing cr, as shown in Fig. [51 
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FIG. 9: The influence of the particle size on the interaction 
potential {$). The graph shows the case of a particle "1" 
(with mass lOMo) interacting with particle "2" (with mass 
50A/o) and particle "3" (with mass 200Mo), respectively, for 
two different values of the ^-potential: C,=ll mV (dashed 
lines), (^=15 mV (full lines). As can be seen, particle "1" 
preferentially aggregates with particle "2" rather than with 
particle "3", to which corresponds a higher potential barrier. 



5. COMPARISON WITH THE EXPERIMENTAL 
DATA 



Eq. 8 provides a connection between the experimen- 
tally measured values of the radius R and the ^-potential 
of the cluster aggregates and the unknown value of the 
variance a"^ of the surface potential V'- This quantity can 
be evaluated for the different systems shown in Fig. [2 i.e., 
cationic particles in the presence of anionic polyions and, 
conversely, anionic particles in the presence of cationic 
polyions. Fig. [TUl shows the dependence of the standard 
deviation a normalized to the value of the ^-potential, 
jcr/CI, as a function of the normalized molar charge ratio 
f/^o- Here, ^ is defined as the ratio between the polyion 
and lipid molar concentrations and ^o is the value of xi 
at which the C-potential goes to zero. As can be seen, 
there is a more or less pronounced decrease of |cr/C| as 
a function of charge ratio ^/^o for all the systems inves- 
tigated, as one can expect when more and more polyion 
adsorption results in a more uniform charge distribution. 

A final comment is in order. Experimentally, in some 
of the colloidal systems investigated, it has been observed 
that, close to the point of charge inversion, the aggregates 
do not reach an equilibrium radius, but their size contin- 
ues to increase with time, until, in long time limit, aggre- 
gates flocculate. This window of ^/^o values is indicated 
as a dashed area in Fig. [TOl Interestingly, if the value 
of the aggregate radius measured just before flocculation 
is used in Eq. [5] in conjunction with the corresponding 
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FIG. 10: The ratio \u/C,\ calculated from eq. 8 on the basis of 
the experimental values of radius 71 and the i^-potential of the 
aggregates formed in different polyion-induced particle aggre- 
gations. The marked regions correspond to the experimental 
instability of the aggregates, where their average size increases 
with time, until, in long time limit, they flocculate. A: Posi- 
tive charged particles: (v)- DOTAP liposomes (0.8 mg/ml) 
and polyacrylate sodium salt; (•): DOTAP liposomes (1.7 
mg/ml) and polyacrylate sodium sall*^. B: Negative charged 
particles: (o): hybrid niosomes and a-polylysine; (A): hybrid 
niosomes and e-polylysine; (□): hybrid niosomes and PEVP 
(ionization degree 65%); (o): hybrid niosomes and PEVP 
(ionization degree 95%)ii. Hybrid niosomes are built up by 
Tween20, Cholesterol and dicethylphosphate and the cationic 
polyion PEVP is Poly[N-ethyl-4-vinyl pyridinium] bromide. 



C-potential value, no physical solution for <j is recovered. 
More precisely, the region of R and C, values for which 
no physical solutions for a exist is delimited by the curve 
a — Q. Along this curve, the equilibrium radius depends 
on ^-potential according to 



R 
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(11) 



. All the points falling in the dashed windows (instability 
windows) of Fig. [10] are positioned in this area oi R— C 
plane (dashed area in Fig. Ilip thus we find particularly 
convincing the fact that the Vclcgol and Thwar potential 
provides an estimate of this unstable region because of 
the close agreement with the experimental results. 



6. CONCLUSIONS 

A large body of experimenta l "'^'^'^° , theoretical^!^ and 
simulation22, investigations have shown that linear flex- 
ible or semiflexible polyions induce the aggregation of 
oppositely charged colloidal particles displaying a large 
variety of possible structures, thus emphasizing the many 
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FIG. 11: Prohibited region under the curve defined by ex- 
pression (|lip for T — 298 K and permittivity e of water. 



facets of charged macroion complex formation. Among 
these structm^es, clusters of liposomes stuck together by 
oppositely charged polyions form a class of model col- 
loids in soft-matter physics which are of the great im- 
portance in many biotechnological implications. In these 
cases, the attractive interaction contribution to the inter- 
particle potential is originated by a correlated adsorption 
at the particle surface, which causes a non-homogeneous 



charge distribution at the particle surface. 

Since the major driving force of these processes is 
of electrostatic origin, more refine interaction poten- 
tials should be used to take into account some pecu- 
liar features in the aggregate cluster formation. We have 
modeled particle interactions by means of a mean force 
inter-particle potential recently proposed by Velegol and 
Thwar, who developed a closed form analytical model to 
estimate the effect of a non-uniform charge distribution. 
Within the scenario of polyion-induced charged parti- 
cle aggregation, this potential provides a justification for 
why an attractive contribution arises in like-charged par- 
ticle interactions in the presence of heterogeneities on the 
charged surface. 

Monte Carlo simulations have been carried out for an 
ensemble of particles interacting via a DLVO-like- Velegol 
and Thwar potential, for different values of the parame- 
ters characterizing the system. Simulations qualitatively 
reproduce all experimentally observable trends in a vari- 
ety of different colloidal systems, from positively charged 
liposomes in the presence of anionic polyions to nega- 
tively charged hybrid niosomes in the presence of cationic 
polyions. In particular, simulations clearly evidence the 
formation of an arrested cluster phase, which is the sign 
of a kinetic arrested state in low-density colloidal suspen- 
sions. 

For electrostatically highly coupled systems, the use 
of the Velegol and Thwar potential offers an interesting 
promise for the understanding of the role of electrostatic 
attractive and repulsive interactions between charged col- 
loidal particles and oppositely charged polyions. 
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